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Abstract. An experimental method to determine the (n,y) cross section of short-living s-proccss 
branching points using data of the inverse (y,n) reaction is presented. The method was used to 
\ observe the branching point nucleus 95 Zr because the elemental abundance patterns corresponding 

to this branching point cannot be reproduced by full stellar models and a possible error source is the 
neutron capture cross section of 95 Zr. The analysis of the experiment is still under progress, we will 



outline the current status in this manuscript. 
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g : 1. INTRODUCTION 

The nucleosynthesis of the elements heavier than iron is today mainly explained by three 
processes: s-, r-, and p-process. s- and r-process are based on neutron-capture reactions 
with adjacent /3 -decays while the p-process is governed by photodisintegration reactions 
such as (y,n), (y,p), or (y,a). 

The distinction between s- and r-process is due to the different neutron densities n n 
involved. During s-process nucleosynthesis n n is in the order of 10 8 neutrons/cm 3 while 
typical r-process sites deal with n n > 10 20 neutrons/cm 3 [1]. Hence, the neutron capture 
rates A.( n ,y) are larger than typical /3 -decay rates during r-process nucleosynthesis 
(\ n ,y) ^ ^*j3> r: rapid). In the s -process the situation is the other way round (^( n ,y) ^ ^j3> 
s: slow) and the involved nuclei are close to the valley of /3 -stability. 

However even in a s-prozess scenario, if the half-life T\n is long enough and the neu- 
tron capture cross section high enough another neutron capture might take place and 
the path "branches" out. Therefore, these nuclides are called branching points of the s- 
process. The branching ratio - i.e. how often each of the paths is taken - determines the 
corresponding elemental abundance patterns. Using a model for ^-process nucleosyn- 
thesis and knowing precisely the nuclear physics input (half-life T\n and Maxwellian 
averaged capture cross section (MACS) at typical s-process temperatures kT = 30 keV) 
it is possible to determine the astrophysical parameters temperature T and neutron den- 
sity n n . 

In the so-called classical approach temperature T and neutron density n n are consid- 
ered to be constant. Thus, three different components produced with different astrophys- 
ical parameters are needed to reproduce the observed abundances of s-only nuclei: the 
weak component corresponds to mass numbers A < 90 while the strong component only 



describes the termination of the s-process path at lead and bismuth. In between these 
borders the isotopes belong to the so-called main s-process [2]. The abundances of s- 
only isotopes that are not affected by a branching can be reproduced with a mean square 
deviation of about 3% using this simple model Jj]. 

A more realistic approach is a full stellar model, e.g. the AGB star model described 
in Ref. H). These models need very precise nuclear physics input data to reproduce the 
5-only abundances. Especially the MACS of several branching points hamper the relia- 
bility of the predictions: due to the lack of experimental data the theoretical predictions 
of the MACS sometimes show a broad spread (see o]). 

A direct measurement of the MACS of the branching points is only possible if their 
half-lives are in the order of years (e.g. 147 Pm with 7\ n = 2.62 yr B6J]). In the case of 
short-living branching points with half-lives of about a dozen days or even less, adequate 
samples are not available. Thus, direct measurements are possible. 

To solve this problem we have investigated an experimental method using the data of 
the inverse (y,n) reaction as described in Sectional Section |3]presents 95 Zr focusing on 
astrophysical aspects as well as on the constraints of our experimental method for this 
nuclide. The present status of the analysis is summarized in Section l3~3l In Section|4|we 
describe the next step in the data evaluation. 



The aim of our experimental method is to constrain the theoretical predictions of the 
MACS of short-living s-process branching nuclei. Short-living under s-process condi- 
tions denotes half-lives of about a dozen days or less, i.e. direct measurements are not 
possible due to the lack of adequate samples. Therefore, our method is based on data of 
the inverse (y,n) reaction. 

The first step is to measure the so-called energy-integrated cross section I a using the 
photoactivation technique: 



with o(E) being the (y,n) cross section and Ny(E,E max ) the experimental photon flux. 

We use the monoenergetic electron beam provided by the S-DALINAC CD to produce 
a continuous bremsstrahlung spectrum by fully stopping the beam in a thick rotating 
copper target. The spectrum ranges from E = to E e , so that E mSLX of Eq. (HJ) equals 
E e . Due to the location of our setup directly behind the injector of the S-DALINAC, 
our maximum energy is limited to about E e - = 10 MeV presently [8]. While the shape 
of N r (E 

? ^max) is well known the absolute value has to be determined by a relative 
measurement with gold as a calibration standard. 

After irradiating the targets the yield Y of the activation is measured by counting the 
y-rays emitted in the decay of the produced radioactive nuclei with a HPGe-detector that 
is well-shielded against background by several layers of lead. The yield Y is proportional 
to the energy-integrated cross section I a with the constant of proportionality being 



2. EXPERIMENTAL METHOD 
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explained in detail in Ref. Jsl]. 



Because the (y,n) cross section of Au is known to behave like 

/ e - S \ °- 545 

a Au (£) = 146.2mb-f-^M (2) 

near threshold [10], a predicted energy dependence of the target's cross section Otarget(£) 
can be normalized to the observed ratio I^ get /I^ u by a correction factor / for each 
measured energy E mSLX . hi doing so, either a parametrization or a theoretical calculation 
of a t arget(£) can be used. If / depends on E max the shape of <7 t arget(-E') is not accurately 
described (compare Eq. [TJ and / can not be used for the further procedure without 
modifications. 

Although the photoactivation technique is very sensitive, several limitations concern- 
ing the activation reaction A X(y,n) A_1 X and the following decay of A_1 X have to be 
taken into account. To test the energy dependece of the cross section Otarget(£) it is nec- 
essary to measure at different energies E max . Because our present setup is limited to 10 
MeV, isotopes with neutron thresholds 5 n > 9 MeV cannot be checked. 

If one wants to use naturally composed targets, the fraction of the isotope of interest 
must not be too small and the activation of other isotopes should not produce too much 
background radiation when the decay of A_1 X is observed. To reach a good peak-to- 
background ratio in the measured spectra only a few y-rays with high branchings should 
be emitted in the decay of A_1 X. Furthermore, the analyzed y-decays should not be part 
of a multi-level decay and their energies Ey should not allow self-absorption in the target 
to avoid additional uncertainties in the analysis. 

If the normalization of a theoretical prediction of o tmge t(E) yields a constant factor 
/, we use / to correct the MACS predicted in the same model with equal parameters. 
This second step of using the correction factor for the inverse reaction is based on the 
detailed balance assumption. In our case, scaling the MACS with a constant factor / 
equals scaling the y-ray strength function. For further discussion see 11 ill . 

Although the corrected MACS is still based on a theoretical prediction, its value is 
constrained by experimental knowledge of the inverse reaction due to the determination 
of /. Thus, it is more reliable than a value purely based on theory. 



3. THE BRANCHING POINT % ZR 

3.1. Position on the ^-process path 

The branching point 95 Zr is located near the borderline between the weak and the main 
component of the s -process Q2Q- Therefore, the prediction of the elemental abundance 
patterns corresponding to this branching is a crucial test for the validity of a full stellar 
model described in ||4|]. 

Several problems concering the prediction of the zirconium abundance patterns mea- 
sured in SiC grains are reported in Ref. II 1211 . The uncertainty in the predicted MACS of 
95 Zr is mentioned as one of the possible error sources due to the wide range of predicted 
values from 23 mb to 126 mb at kT = 30 keV (see fill, ifll. Pi. and Ml Another 




FIGURE 1. The ^-process path around the branching point Zr. The main flow of the ^-process is 
marked with thick lines, thinner lines indicate branches. Grey shaded boxes illustrate stable isotopes with 
the numbers at the bottom of the boxes being their natural abundances. The black box of 96 Mo stresses 
that we have an s-only nucleus. The white boxes display j3 -instable nuclei: dashed boundaries are used 
for /3~-decay, dotted lines are used to indicate j3 + - or e-instable isotopes. In this cases the numbers at the 
bottom of the boxes are the half-lives of the ground-states. 



reason for the deviation between the measured and predicted abundance patterns might 
be the uncertainty in the half-life Tmq of 95 Zr at s-process temperatures. However, as 
confirmed by a recent measurement [ 17] the first excited level is at E = 954 keV and 
thus, is not significantly populated at kT = 30 keV. Hence, the half-life of 95 Zr does not 
depend on temperature under s-process conditions and can be omitted as an error source 
because of its small error: 7/ 1/2 ( 95 Zr) = (64.032 ±0.006) d d. 

3.2. Constraints of the experiment 

To use our experimental method the constraints described in Section El have to be 
fullfilled for the activation reaction 96 Zr(y,n) 95 Zr. The threshold of this reaction is 
S n = 7854 keV, which can be easily reached at our setup providing energies up to 10 
MeV. The abundance of 96 Zr in a naturally composed target is 2.8% and thus, high 
enough for the activation method. 

The properties of the decay of 95 Zr are summarized in Fig. |2j The characteristics 
claimed in Section[2]are achieved for Ey — 724.2 keV and Ey = 756.7 keV. 

Due to the half-life of 95 Zr of about 64 days it is mandatory to observe the decay for 
several hours to reach good statistics. During this time the /3 -decay of the daughter 
nucleus 95 Nb also starts to take place. The assignment of the measured y-rays to the 
decays of 95 Zr or 95 Nb can either be done by a determination of the observed half-life 
or by comparing their different behaviour in time. As shown in Fig. |3] the ratio of the 
two peaks at 724.2 and 756.7 keV remains unchanged if the measurement takes place 
different times after activation. Thus, both lines correspond to the same decay and are 




FIGURE 2. Section of the decay scheme of 95 Zr based on 1 18] 



identified to be the most prominent lines of the decay of 95 Zr. In contrast, the peak at 
765.8 keV becomes more and more dominant with time. This behaviour is consistent 
with the one expected of a line corresponding to the decay of the daughter nucleus 
95 Nb. Actually, the energy E r = 765.8 keV is exactly the energy supposed for the most 
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FIGURE 3. Typical decay spectra of 95 Zr each measured for t — 1 8 h. From the top to the bottom panel 
the time between activation and start of the measurement is changing as indicated on the top right. For 
further explanation see text. 



prominent line of this decay. 



3.3. Current status of analysis 



The photoactivation of zirconium has been done at six different energies E max ranging 
from 8325 to 9900 keV. The lower limit is due to the separation energy of the used 
calibration standard: 5 n (Au) = 8071 keV. In a first approach, we used a parametrization 
that describes the energy dependence of the (y,n) cross section near threshold: 



o(E) = o 



E-Sn 



(3) 



In this case, the exponent p changes the shape of the energy dependence, whereas the 
factor Oo is responsible for the absolute value of o(E). Therefore, Go corresponds to the 
correction factor / of Section El Due to the fact that this parametrization is derived from 
the time reversal symmetry of (n,y) and (y,n) reactions, the exponent p is connected to 
the wave-character of the emitted neutron: an exponent of p = 0.5 + / corresponds to an 
Z-wave decay. Thus, exponents less than 0.5 cannot be explained in this basic approach. 

However, if a least-mean-square-fit of our data to this parametrization is carried out 
using Go and p as free parameters, we obtain p = 0.339 and Gq = 31.81 mb as shown 
in Fig.|H The single results <Jo(E mSLX ) fluctuate around the mean value, thus, the overall 
energy dependence of o(E) seems to be described accurately. 
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FIGURE 4. Result of a least-mean-square-fit of Eq. Q to the measured data at different energies E max . 
The dots are the results for Ob derived with an exponent of p = 0.339 while the dashed line is the mean 
value Ob = 3 1 .8 1 mb. The horizontal error bars are due to uncertainties in the beam energy. 

The result p < 0.5 can be explained by assuming a resonance on top of the smooth 
behaviour described by Eq. ©. We measure the energy integrated cross section I a and 
compare it with the result of Eq. ([TJ using o(E) of Eq. ©. Thus, we average the 
measured (y,n) cross section <j(E) over the energy interval ranging from 5 n to E max . 
A resonance in this energy region would yield the measured bigger I a . 

The smaller p, the steeper is the rise of the cross section at the threshold energy S n . 
Hence, the calculated I a becomes bigger and suits the measured value. Thus, we accept 
the derived values for p and Go as a description of the mean energy dependence of the 
(y,n) cross section: g(E) = 31.81 mb ■ ((£ -5 n )/5 n ) a339 . 

Fig. |5] shows the mean g(E) and a theoretical prediction of the (y,n) cross section 
Gtheo(E) based on Hauser-Feshbach calculations Ill9l. l201.lf the shapes of the two energy 
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FIGURE 5. Comparison of the energy dependence of the cross section o(E) for the reaction 96 Zr(y,n). 
The solid line is a theoretical prediction based on Hauser-Feshbach calculations. The dashed line is the 
result of a least-mean-square-fit of the parametrization described in Eq. to our data. 



dependences are compared, a significant deviation is obvious: the theoretical prediction 
rises very steeply close to the threshold, but its slope becomes flatter 50 keV above 
threshold. In contrast, our parametrization is not so steep near threshold but more 
continuously rising. 

Although the parametrization of o{E) describes only a mean energy dependence 
of the cross section, the different behaviour illustrates that the energy dependence of 
the theoretical prediction does not suit the data as well. Obviously, no resonances are 
predicted by theory. A correction factor /, derived as explained in Section |3 is only a 
mean correction. / can compensate the absence of resonances in theory. Furthermore, 
/ is a function of E m?lX . Therefore, the MACS of 95 Zr predicted by this theory cannot 
simply be corrected by this mean value of / but the (n,y) cross section must be corrected 
by f{Ey) before calculating the MACS. f(Ey) can be calculated weighting f(E max ) by 
o(E) - Ny{E,E max ). For this purpose, it is necessary to measure the (y,n) cross section 
as close as possible to the reaction threshold to locate the position of the resonances. 

So far, the experiment was limited by the threshold of the calibration standard. The 
lowest measured energy was E max = 8325 keV i.e. still about 500 keV above the 
threshold 5 n of 96 Zr. In order to measure as close as possible near the target's threshold 
the neutron separation energy of a new calibration standard has to be lower. 

In addition to a low neutron separation energy the calibration standard has to fullfiU 
other characteristics: despite the general features mentioned in Section El the half-life 
Tj j 2 of the produced radioactive isotope has to be short so that the counting can be done 
in short times with good statistics. 

We decided to use 187 Re with a neutron separation energy of only 5 n = 7359 keV and 
a natural abundance of 62.6%. The half-life of the produced 186 Re is 7\ n = 3.72 d and 
allows short measuring times although the branching of the two dominant y-decays is 
quite low: E y = 137.2 keV with 9.42% and E y = 122.6 keV with 0.6%. Self-absorption 
in the rhenium target limits the amount of material, however, based on experience of 
former experiments 111 ill very thin rhenium targets with a thickness of d = 50 /im should 
provide enough activity after a typical activation. 




E^/keV 



FIGURE 6. Comparison of the decay spectra of the calibration standards gold and rhenium. The inlays 
show the 7 peaks that are used for the calibration. The peak-to-background ratio is in the same order of 
magnitude for both standards. 



187 Re was measured in May 2003 and the analysis is still under progress. Fig. |6] 
compares the spectra of rhenium and gold measured after an activation at E max = 
8325 keV. The upper panel shows the spectrum of rhenium measured for three hours 
directly after the activation. The peaks between 750 and 1000 keV correspond to the 
decay of 184 Re that is produced by 185 Re(y,n). Due to these lines the two peaks of 
the decay of 186 Re at Ey = 137.2 keV and Ey = 122.6 keV are located above a high 
background. However, as shown in the inlay the peak-to-background ratio of both lines 
is good enough to analyze the peak volumina with a low statistical error. For comparison 
the inlay in the lower panel shows the two most prominent lines of a gold spectrum 
also measured for three hours. The peak-to-background ratio is in the same order of 
magnitude. 



4. SUMMARY AND OUTLOOK 

The precise knowledge of the (n,y) cross section of the ^-process branching point 95 Zr 
is a crucial input parameter to test the reliability of the AGB star model. Due to the very 
wide range of the theoretically predicted values a measurement observing the inverse 
reaction with the photoactivation technique was carried out at the S-DALINAC in June 
2002. 

During the analysis of the data several problems occured: The energy dependence 



of the (y,n) cross section of 9 °Zr extracted from the data differs significantly from the 
theoretical prediction. In our opinion, the deviation is due to resonances that are not 
included in the used theory. Therefore, the experimental method described in Section E] 
determines an energy dependent correction factor f(E r ) that cannot be used to correct 
the predicted MACS of 95 Zr without further studies. 

Hence, the experiment was extended in May 2003 by measuring the (y,n) cross section 
of 96 Zr closer to the reaction threshold. Therefore, a new calibration standard with lower 
neutron separation energy was measured during the same beam-time. The analysis of the 
new standard 187 Re is still under progress, preliminary results like the spectrum shown 
in Fig.|6]are encouraging that it will work as well as our first standard 197 Au. 

Once a description of the (y,n) cross section of 187 Re near threshold is available 
the new data point of 96 Zr taken at E max = 8100 keV has to be added to the existing 
analysis. Therewith, the discrepancy between the experimental result and the theoretical 
prediction will hopefully disappear or at least become explainable. The outcome of the 
whole analysis that is the experimentally confirmed MACS of 95 Zr can then be used in 
the AGB star model and might solve the problems in explaining the observed zirconium 
abundance patterns in SiC grains. 

In future, we will take more data points at different energies E mSLX in order to reduce 
difficulties due to resonances on top of the cross section. If the distance between two 
data points AE max is smaller, the determination of the location of a resonance is possible 
with higher precision. 
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